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We present experiments on the driven dynamics of a two-level superconducting artificial atom. 

The driving strength reaches 4.78 GHz, significantly exceeding the transition frequency of 2.288 GHz. 

The observed dynamics is described in terms of quasienergies and quasienergy states, in agreement 
with Floquet theory. In addition, we observe the role of pulse shaping in the dynamics, as determined 
by nonadiabatic transitions between Floquet states, and we implement subnanosecond single-qubit 
operations. These results pave the way to quantum control using strong driving with applications 
in quantum technologies. 


Monochromatic driving is the most common tool in 
quantum control, applicable to various physical systems 
including nuclear and electronic spins, atoms, ions, su¬ 
perconducting qubits, and quantum dots [1]. For driv¬ 
ing that is weak compared to the relevant transition fre¬ 
quency, the dynamics can be described in terms of Rabi 
oscillations between energy eigenstates. In contrast, with 
strong driving the commonly used rotating wave approx¬ 
imation [2] breaks down, resulting in complex evolution. 
Strong driving dynamics is most adequately described 
in the framework of Floquet theory [3], where the state 
of a driven system is expressed in terms of quasiener¬ 
gies and quasienergy states. Exploring this more general 
framework expands the field of quantum control, gaining 
increasing relevance as current experiments on the imple¬ 
mentation of high-fidelity quantum gates [4] and protec¬ 
tion against decoherence [5] are performed with a driv¬ 
ing strength that is a significant fraction of the transition 
frequency. In addition, strong driving is relevant in the 
fields of quantum sensing, for phase measurements [6], 
and quantum simulation, for designing effective Hamil¬ 
tonians in the emerging field of Floquet engineering [7]. 

In this Letter we report experiments on the dynam¬ 
ics of an artificial atom, a superconducting quantum 
bit [8, 9], strongly driven by a microwave field. Strong 
driving has been studied in the field of atomic physics, 
using either optical [10] or radio-frequency pulses [6]. In 
experiments with NV centers in diamond, time dynam¬ 
ics was observed for driving strength up to values com¬ 
parable to the transition frequency [11]. Superconduct¬ 
ing qubits display a naturally strong coupling to electro¬ 
magnetic fields due to their mesoscopic character. Pre¬ 
vious experiments on strong driving of superconducting 
qubits have mostly addressed the steady-state response 
to continuous waves [12-22]. A few experiments have ob¬ 
served time-domain Rabi oscillations [23-25] with a driv¬ 
ing strength exceeding the transition frequency, and two 
of these demonstrated good agreement with the theoret¬ 
ically predicted Bessel-function dependence of the Rabi 


frequency [23, 25]. In our experiments, we use quantum 
state tomography to investigate the dynamics of a su¬ 
perconducting qubit strongly driven by microwave pulses 
with controllable shape. The observed system dynam¬ 
ics is very well described in terms of quasienergies and 
quasienergy states, as predicted by Floquet theory. In 
particular we observe several frequency components in 
the dynamics, in very good agreement with theory. We 
find that the switching on and off of the driving pulse 
plays an important role in the qubit evolution, as de¬ 
termined by adiabaticity conditions in the Floquet pic¬ 
ture [26]. We also used strong driving for fast, sub¬ 
nanosecond, preparation of qubit states. 

The artificial atom in our experiment is a supercon¬ 
ducting flux qubit [27]. Among the different types of su¬ 
perconducting qubits, flux qubits have the advantage of 
high-level anharmonicity, leading to ideal two-level sys¬ 
tem behavior, and of strong coupling to electromagnetic 
fields [28, 29] which enable strong driving. Qubit state 
measurement is performed by probing microwave trans¬ 
mission through a resonator coupled to the qubit (see 
Fig. 1(a)), in the dispersive regime of circuit quantum 
electrodynamics [30]. 

The qubit Hamiltonian is given by H{t) = —^Uz — 
in a basis formed by symmetric and antisymmet¬ 
ric combinations of clockwise and anticlockwise persis¬ 
tent current states in the qubit loop [27]. Here A, the 
minimum energy level splitting, is a fixed parameter, and 
e{t) = — 4)o/2), with the magnetic flux ap¬ 

plied to the loop dependent on the time t. The magnetic 
flux $(t) = <I>s + 4>d(t) with a static flux generated 
by a superconducting coil and ^d{t) a time-varying mag¬ 
netic flux coupled to the qubit through a waveguide ter¬ 
minated by an antenna (see Fig. 1(a)). The usual ap¬ 
proach employed to generate control pulses is based on 
using a modulator to shape the quadratures of a contin¬ 
uous wave produced by a frequency synthesizer. Here 
we use instead a new generation of high-speed arbitrary 
waveform generator (AWG) to directly synthesize the mi- 
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FIG. 1. (a) Schematic representation of the experimental 

setup. The qubit, formed by a superconducting loop in¬ 
terrupted by Josephson junctions (cross symbols), is cou¬ 
pled to a superconducting coplanar waveguide resonator. 
Readout is based on the transmission of a microwave pulse 
from the resonator input (left) to its output port (right). 
A waveguide (bottom) is used to couple microwave control 
pulses to the qubit, (b) Representation of qubit control 
pulses, with rise and fall times G and t{ respectively, and 
maximum-amplitude duration tp. The thick line indicates 
the pulse envelope, which reaches a maximum Am- During 
rise and fall, the envelope is shaped as (1 — cos 
and (1 -I- cos (7r(t — tp — tr)/tt)). (c) Qubit transition fre¬ 
quency tuoi from spectroscopy measurements versus the static 
magnetic flux $ 3 . The continuous line is a fit of the transition 
frequency, yielding the parameters A = 27r x 2.288 GHz and 
Ip = 690 nA. 


crowave pulses [31], leading to the time accuracy required 
for control with subnanosecond resolution pulses. 

A plot of the qubit transition frequency versus the 
static flux obtained by spectroscopy with weak and 
long microwave pulses, is shown in Fig. 1(c). All the ex¬ 
periments reported in this Letter are performed at the 
symmetry point ($s = ^“ 0 / 2 ), where the qubit transition 
frequency wqi = A = 27 r x 2.288 GHz. We use amplitude 
shaped pulses e{t) = 2A(t) cos(wt), with A{t) character¬ 
ized by a maximum amplitude Am and rise and fall times 
denoted by tr and tf respectively (see Fig. 1(b)). At the 
symmetry point, the energy relaxation and pure dephas¬ 
ing times are given by Ti = 1.8 /iS and iRamsey = 0.3 /xs. 
These coherence times, currently limited by quasiparti¬ 
cle tunneling, microscopic two-level systems, and charge 
noise [32, 33], can be further improved by infrared shield¬ 
ing techniques and improved qubit design without im¬ 
pairing the ability to strongly drive the qubit. 

Experiments are performed by repeating, typically 
16,384 times, a sequence formed of state reset, control 
using an applied pulse, and measurement in the energy 
eigenbasis. Figure 2(a) shows the qubit’s average excited 
state probability versus the duration of the microwave 
pulse, with driving on resonance. The waveform is de¬ 
fined with zero rise and fall times; however, the actual rise 


and fall times are determined by the analog bandwidth of 
the AWG and are specified to be shorter than 22 ps [31]. 
For weak driving (Fig. 2(a), top panel), sinusoidal oscil¬ 
lations are obtained, as predicted based on the rotating 
wave approximation. With a large Rabi driving strength 
(Fig. 2(a), bottom panel), large amplitude oscillations are 
accompanied by smaller amplitude faster oscillations, a 
signature of non-negligible counter-rotating term effects. 
The different frequency components are clearly visible in 
the Fourier transform of the signal (Fig. 2(b)). For a wide 
range of the driving strength from 27rx0.20 GHz to 
27rx4.78 GHz, the Fourier transformed data are shown 
in Fig. 2(c). 

The presence of the various frequency components in 
the Rabi oscillations can be understood based on Floquet 
theory, which predicts that for a time-periodic Hamil¬ 
tonian with period T the quantum state is given by 
IV’(^)) = 1 with Cj the quasienergies 

and \uj(t)) the quasienergy states, periodic in time with 
period T. As a result, the probability to find the sys¬ 
tem in its excited state is expected to show oscillatory 
behavior with frequency components nco and ±Ae -I- nw, 
with Ae the quasienergy difference, uj = 27r/T the driv¬ 
ing frequency, and n any integer number. The harmonic 
drive signal used in our experiment has the additional 
symmetry e{t + T/2) = —e{t), and as a result only com¬ 
ponents with even n values are present [34]. Fig. 2(d) 
shows the extracted frequency components versus driv¬ 
ing amplitude. We compare the experimental results 
with calculations of the quasienergies based on numer¬ 
ical simulations (solid lines) and an analytical expres¬ 
sion (dashed lines). The latter, obtained based on ap¬ 
proximate diagonalization after transformation to a ro¬ 
tating frame [31], gives a quasienergy difference Ae = 

uj\J (1 — Jo (^))^ + Ji with Jg/i Bessel functions 
of the first kind and order 0/1. This formula provides 
a good approximation for the case of a two-level system 
biased at its symmetry point and driven on or near res¬ 
onance with arbitrary strength, complementing previous 
theoretical work where the weak- and strong-driving lim¬ 
its of this formula had been derived [21, 22]. Additional 
experiments were performed with the qubit driven off- 
resonance, with a driving frequency w = 27r x 1.373 GHz 
(see Fig. 2(e) and (f)). The Fourier transform of the qubit 
population signal and the identified frequency compo¬ 
nents are shown in Figs. 2(g) and (h) respectively. Good 
agreement with the predictions of numerical calculations 
is observed in this case as well (see Fig. 2(h)). 

Tomography experiments confirm the role of the 
counter-rotating terms in the driven evolution of the 
qubit. Fig. 3 shows results of state tomography versus 
the duration of driving pulses for two values of the driving 
amplitude. Am = 27rx0.10 GHz and Am = 27rx0.46 GHz, 
and zero rise and fall times. For both values of the 
driving amplitude, high-amplitude oscillations are ob- 
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FIG. 2. Coherent oscillations versus driving amplitude for resonant driving (uj — A, (a-d)) and off-resonance driving (oj = 
0.6 X A, (e-h)). (a,e) Qubit excited state probability Pi versus control pulse duration tp for Am ~ 0.10 and 1.00 GHz (a) 
and Am = 0.30 and 1.44 GHz (e). (b,f) Discrete Fourier transforms of the oscillations in (a) and (e) respectively. (c,g) Color 
plot of the Fourier transform of population oscillations versus frequency and driving pulse amplitude. (d,h) Positions of peaks 
in the Fourier transform of coherent oscillations versus driving amplitude extracted from the data in (c) and (g) (dots). The 
lines are plots of nui, nui — Ae, and nuj -F Ae respectively, with the quasienergy difference Ae determined numerically and n 
an even integer. For the resonant driving case (d) we show in addition corresponding curves (dashed) based on the analytical 
approximation for Ae discussed in the text (see also Ref. [31]). 


served, ■with a period corresponding to the quasienergy 
difference. In the weak-driving limit, these oscillations 
are the usual Rabi oscillations. High-frequency compo¬ 
nents are observed in addition, with a significant am¬ 
plitude at strong driving, reflecting the presence of the 
non-negligible counter-rotating wave component. The re¬ 
sults of tomography experiments are in very good agree¬ 
ment with predictions of numerical simulations of the 
Schrodinger equation (see Fig. 3). 

The presence of the fast oscillatory terms in the 
driven evolution depends not only on the pulse ampli¬ 
tude, but also on the pulse turn-on and turn-off times. 
Fig. 4(a) shows qubit state oscillations for a driving 
strength = 2tt x 1.33 GHz, and different rise and fall 
times. Fast oscillatory terms are gradually suppressed 
as the turn-on and turn-off times are increased. We 
emphasize that fast oscillatory components in the oscil¬ 
lations are completely suppressed for slow pulse turn¬ 
on and turn-off despite the fact that during most of 
the driven evolution the driving amplitude is compa¬ 
rable with the transition frequency. The absence of 
fast oscillations for slow turn-on and turn-off can be 
understood based on adiabaticity in the Floquet pic¬ 
ture [26]. Indeed, the time-dependent qubit state can be 
written, up to an overall phase and a geometric phase, 
as \'ip{t)) = co(t)|uo(x4,t))-F ci(t)e“*-^o 
with |uo(H., t)) and |Mi(H,f)) the instantaneous driving- 
amplitude-dependent quasienergy states. The initial val¬ 
ues of the coefficients Cq and ci are determined by the rep¬ 
resentation of the initial qubit state, which is the ground 
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FIG. 3. Measurements and simulations of the evolution of 
the Bloch vector components, given by the average values 
of the Pauli Ua (a = x, y, x) operators, after a pulse with 
zero rise and fall time, versus the length of the pulse for 
Am = 27rx0.10 GHz (a) and Am = 27rx0.46 GHz (b). The ex¬ 
perimental results are in excellent agreement with the results 
of numerical simulations. 


state, in the basis formed by the states |mo,i(0,0)) = 
(jo) ± |l))/-\/2, with |0)(|1)) the ground/excited state of 
the qubit [31]. For slowly varying driving amplitude A{t)^ 
the evolution is adiabatic in the Floquet basis, and there¬ 
fore the coefficients cq and Ci maintain their initial values. 
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The dynamics of the qubit in this case, using the Bloch 
sphere representation (see Fig. 4(d)), can be understood 
by the rotation of the pseudospin representing the state 
around a fictitious field determined at any given time by 
the difference Ae between the quasienergies. Similarly to 
the weak-driving case, the qubit simply undergoes a Rabi 
rotation between the ground and excited states, with a 
rotation angle given by 1/2 Ae{t)dt. 

With short rise and fall times, the evolution of the 
qubit at the beginning and end of the pulse is nonadia- 
batic in the Floquet representation. Nonadiabatic effects 
can be described by unitary transformations Up,rise and 
UpMi the beginning and the end of the pulse respec¬ 
tively (see Fig. 4(e)). The latter depends periodically on 
the pulse duration with period T, leading to fast oscilla¬ 
tions of the final state of the qnbit. For a driving ampli¬ 
tude Am = 27rx 1.33 GHz, the qubit population dynamics 
is well described by a sum of oscillatory terms at frequen¬ 
cies Ae, 2oj -I- Ae, and 2oj — Ae. In Fig. 4(c) we plot the 
amplitude of the high-frequency components, at 2a;± Ae, 
versus the pulse rise and fall time. The experimental re¬ 
sults are in good agreement with values extracted based 
on numerical simulations of the qubit evolution. 

In additional experiments (see Fig. 4(b)), we observed 
the evolution of the qubit with strong pulses and asym¬ 
metric rise and fall times. The final state of the qubit 
displays fast oscillations for pulses with slow rise and fast 
fall, whereas fast oscillations are absent for pulses with 
fast rise and slow fall. This observation confirms the 
asymmetric role of the two rotations, Upyise and llpjaii- 

We next discuss the use of strong driving for fast quan¬ 
tum gates, specifically qubit state preparation. Starting 
with the qubit in its ground state, we apply pulses with 
a driving strength Am = 27r x 0.46 GHz and rise and fall 
times of approximately 20 ps, defined by the AWG band¬ 
width. The state (|0) — z|l)) /V2 is prepared in 0.48 ns 
with a fidelity of 0.9996 ±0.0006 [31]. Similarly, state |1) 
is prepared in 1.08 ns with a fidelity of 0.9969 ± 0.0008. 
We have performed numerical simulations of state evolu¬ 
tion, which predict state preparation fidelities of 0.9997 
and 0.9976 for states (|0) — zjl)) /V^ and |1) respectively, 
in good agreement with the experimental results. Future 
work should address the optimization of gate fidelities 
based on randomized benchmarking [4, 35, 36]. 

Our work demonstrates the feasibility of using strong 
driving for the control of superconducting artificial 
atoms. The dynamics was analyzed in the framework of 
Floquet theory. The consideration of adiabaticity in the 
Floquet picture provides a valuable viewpoint on dynam¬ 
ics, applicable well beyond the regime where the rotating 
wave approximation holds. Our experimental demon¬ 
stration brings very exciting prospects for experiments 
addressing the interplay between Floquet dynamics and 
environmental effects [37, 38]. We expect that our results 
will stimulate new work across a broad range of fields 
including quantum computing, open system dynamics. 


quantum simulation, and quantum sensing. 
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FIG. 4. (a,b) Qubit excited state probability Pi versus pulse duration tp for various rise (fall) times tr (tf) and equal maximum 
amplitude Am = 2n x 1.33 GHz. Panels (a) and (b) show data with symmetric or asymmetric rise and fall, (c) The measured 
(dots and squares) and simulated (continuous and dashed lines) fast oscillation amplitudes, at frequencies 2tu — Ae and 2uj + Ae 
respectively. (d,e) Adiabatic (d) and nonadiabatic (e) evolution in the Floquet picture. The state is represented on the Bloch 
sphere, with a basis chosen such that the instantaneous quasienergy states |uo,i) are in the equatorial plane. The state vector 
evolution (thin arrows) is a rotation around a fictitious field (thick arrows). The initial qubit state is (|mo) + |wi))/\/2. In 
the adiabatic case (d), the state evolution is described by a phase — fg Ae(t)dt applied to |ui), which is equivalent to rotation 
around the fictitious field Ae{t) aligned with |uo). In this picture, the evolution of the qubit is a simple rotation, although in 
the energy eigenbasis the qubit undergoes complex dynamics. In the nonadiabatic case (e), transitions between Floquet states 
arise during pulse turn-on and turn-off, characterized by the unitary transformations Uf, rise and [/F,faii respectively, which 
correspond in general to rotations around axes that are not parallel to Ae{t). The asterisk {*) indicates that the actual rise or 
fall time is around 20 ps, determined by the analog bandwidth of the AWG. 
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Supplemental Material 

EXPERIMENTAL METHODS 
Sample fabrication and parameters 

The device (see Fig. SI for images) is fabricated on a high-resistivity silicon substrate, in two steps. Firstly, the 
resonator and the control lines are defined by optical lithography, followed by evaporation of a 190 nm thick aluminum 
layer and liftoff. In the second step, a bilayer resist is patterned by electron-beam lithography. After an argon milling 
step, shadow evaporation of two aluminum layers, 40 and 65 nm thick respectively, followed by liftoff, defines the 
qubit. The resonator is formed by a coplanar waveguide interrupted by coupling capacitors at the input and output 
ports. The fundamental half-wavelength mode used for qubit readout has a resonance frequency Wr = 27r x 6.641 GHz. 
The qubit and resonator are strongly coupled using the inductance of a shared line, with a coupling strength for the 
fundamental mode g = 2 tt x 537 MHz [SI, S2]. 


Measurement setup 

The experiments are performed in a dilution refrigerator at a temperature of 35 mK. The fabricated device is 
enclosed in a copper box, which is placed inside a three-layer high-permeability metal shield. An active magnetic 
field compensation system is used to further reduce the effect of fluctuations of the external magnetic field. The 
device is connected to room temperature electronics using coaxial cables, which include attenuators and filters placed 
at different temperature stages. The signal at the output port of the resonator is amplified using a low-noise high 
electron mobility transistor amplifier with a noise temperature of 4 K. The qubit control pulses are directly synthesized 
by an arbitrary waveform generator Tektronix AWG70002A at a sampling rate of 25 GS/s. The AWG has an analog 
bandwidth of 13.5 GHz and an intrinsic rise/fall time of less than 22 ps. 


Pulse calibration for quantum state tomography 

The quantum state tomography experiments are performed by applying pre-rotation pulses to the qubit before the 
readout. We use pre-rotation pulses denoted by unitary operations including identity / and 7r/2 rotations around the 
X- and j/-axis, denoted by Rx{Tr/2) and Ry{n/2) respectively. The pre-rotation pulses Rx{Tr/2) and Ry{n/2) are defined 
with tr = tf = 0.2 ns and Am = 27r x 130 MHz, where the dynamics are well described by the RWA. The length and 
the phase of these pulses are chosen so that their rotation angles, as well as the angle between their rotation axes, are 
calibrated to tt/2. To calibrate the rotation angles and the rotation axes of these pulses, we manipulate the qubit with 
a pulse sequence consisting of the pre-rotation pulses and then measure the final state of the qubit. For calibrating 
the rotation angle 0 of pulse Rx(y){0), we use a pulse sequence corresponding to a unitary operation (y) 
which amplifies the rotation angle error by 2n -I- 1 times and projects this angle error to the measurement basis. 



FIG. SI. (a) Scanning electron microscope (SEM) image showing a qubit embedded in the resonator. An antenna (right) is 
used to couple microwave control pulses to the qubit, (b) SEM image of a qubit device nominally identical to that used in this 
work. 
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Using the above pulse sequence with n = 5, we determine an optimized pulse length of tp = 1.9 ns and an upper 
bound of 0.003 rad on the error of the rotation angles around both axes. For calibrating the rotation axes, we use a 
pulse sequence corresponding to unitary operation Ryi (f) 1 [Rx (f)] [Ry' (f)] [ Rx{^), where Ry/ (f) denotes 
the rotation around the axis y' to be calibrated. Using pulse sequence with n = 5 which amplifies the error in the 
rotation axis by 2n times, we determine an upper bound of 0.002 rad on the error of the rotation axis y'. These 
rotation-angle and rotation-axis errors lead to an error less than 0.003 for each Bloch vector component (7^. These 
errors are much smaller than the statistical errors from the measurements which will be discussed in the next section. 


Quantum state preparation with strong pulses 

The density matrices of states are reconstructed from the quantum state tomography data using maximum likelihood 
estimation. A direct comparison between the reconstructed state p from experimental results and the ideal state pideai 
is given by the state fidelity F = Tr (i/VPi^gai/^V^ideai) ’ Statistical errors of the reconstructed states and their fidelities 
are determined by the parametric bootstrapping method [S3]. The procedure is described as follows: 1. Calculations 
of the standard deviation a of the Bloch vector a = (ax, ay,az) obtained from the tomography experiments using the 
measurement statistics. According to the central limit theorem, the average of the random variables ai = ai^n/n 
(i = x,y,z) has a normal distribution with a mean of (ui) and a standard deviation of a(ai)l^pn with n the number 
of repeated experiments. 2. Generation of random data of the Bloch vector from the normal distributions above. The 
size B of the random data need to be large, usually B > 100. 3. For each random data, the maximum likelihood 
estimation is used to estimate the quantum state p and the state fidelity F with respect to the target state is calculated. 
Finally, statistics for all the data of size B is collected. The variance of data serves as an estimator of the standard 
error. 

Using pulses with = 2tt x 0.46 GHz and rise and fall times about 20 ps, we prepare states targeting |—U) = 
(jo) — i\l))/'j2 and |1) from the ground state |0) with a total pulse length of 0.48 ns and 1.08 ns respectively. The 
reconstructed density matrices of the actual states are determined to be: 


P\-Y) = 

( 0.511048 

-0.0145217-b 0.499667Z 

1^-0.0145217- 0.4996671 

0.488952 

P\l) = ( 

0.00590452 

-0.0709229 -b 0.0289758i 

-0.0709229 - 0.02897581 

0.994095 


of which the fidelities are 0.9996 ± 0.0006 and 0.9969 ± 0.0008 respectively. 


CALCULATION OF QUASIENERGIES AND QUASIENERGY STATES 
Matrix form of the Floquet Hamiltonian 

Here, we derive a matrix form for the Floquet Hamiltonian of a qubit under harmonic driving for numerical 
calculation of the quasienergies. In the derivation below, we shall generally follow similar steps to those given in 
Ref. [S4]. 

We consider a two-level artificial atom under harmonic driving. The Hamiltonian is given by: 

H =—^az + Acos{ijjt) ax- (SI) 

We consider units with h = \ from here onwards. In order to simplify the appearance of expressions that will appear 
in the derivation below, we make a basis transformation described by the operation 

^rot = (S2) 

i.e. a 7r/2 rotation about the y axis. The Hamiltonian is then transformed into the form 

Hrot =- ^ax - Acos(ujt)az- 


(S3) 
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Since we are dealing with a two-dimensional Hilbert space, the time-periodic Hamiltonian gives rise to two Floquet 
states, which can be expressed as 


\'^F,j{t)}=e (S4) 

where Cj are the quasienergies and \uj{t)) are the periodic components in the Floquet states, i.e. time-dependent 
quantum states with period T = 27r/a;: 


bW)= E 




n— — oo 


Substituting the above expression in the Schrodinger equation, we find the relation 


A 


A 


+ nu „) - (|Uj>-l) + \Uj,n+l)) 


The above set of equations can be expressed as a single equation: 

e,\U,) = Hp\U,) 


(S5) 


(S6) 


(S7) 


given by 


5 1,4-5 

n,t 5 


Uj,n+l,l, ■ 
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(n — l)w 
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A 
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A 

(n — l)w 


A 
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0 

0 

0 


1 

0 

mo 

^A 
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2 

2 



0 

A 

A 
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nuj 


0 

0 

A 

9 

0 

(n -1- l)a; 
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A 
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0 

0 

0 
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0 

0 

0 

A 

\ 

2 

(n -I- l)w 


(S8) 
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One can obtain a good approximation for the quasienergies and quasienergy states by numerical diagonalization of a 
truncated version oi Hp as long as the truncated matrix remains sufficiently large. In our calculations, we truncate 
the matrix Hp with n ranging from —50 to 50. 


Choice of quasienergies and quasienergy states 

Although the eigenvalue problem in Eq. (S7) has an infinite number of solutions, there are only two inequivalent 
solutions and all the other solutions are copies obtained by shifting an integer number of quanta noj between the 
quasienergy and the periodic part of the Floquet states. We take the two eigenvalues of which difference corresponds 
to the Rabi frequency Or = |A — wj in the A —>■ 0 limit as the two inequivalent solutions. The consideration of 
evolution in the Floquet picture requires the decomposition of qubit states into quasienergy states at A = 0 and t = 0. 
While any two orthogonal states are a proper choice of Floquet states at A = 0, we chose |uo,i) = (|0) ± |l))/-\/2. 
This choice corresponds to assuming a finite value of A and then taking the limit A —>■ 0. 


Analytical formula for quasienergies 

In this section, we derive an approximate formula for quasienergies of a qubit biased at the symmetry point and 
subject to harmonic driving. Starting from Eq. S8 in the last section, we now perform a basis transformation that 
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physically corresponds to going to a rotating frame (with a time-dependent rotation frequency). The basis states after 
the transformation are related to those before the transformation by the formula [S5] 

^ )0, To ^ ,0, Ji ^ —^ , 0,... }, 

|u,0,J_i ,0,Jo ,0,Ji (S9) 

The Hamiltonian for the new basis now reads 



(sio) 

It is helpful at this point to use the picture of perturbation theory and think of all the off-diagonal matrix elements 
as a perturbation. Since there are a large number of matrix elements in the perturbation component and many of 
these will have negligible effects, we will try to identify the elements that need to be kept for an accurate description 
of the system. The unperturbed Hamiltonian is given by 



(n-l)w 0 0 0 0 0 

0 (n-l)u; 0 0 0 0 

0 0 nw 0 0 0 

0 0 0 nw 0 0 

0 0 0 0 (n-kl)w 0 

0 0 0 0 0 {n + l)oj 


(Sll) 


Let us now calculate the quasienergy and Floquet state starting from the third entry in the matrix above, i.e. the 
one that corresponds to the first of two appearances of the unperturbed energy nu. In the unperturbed Hamiltonian, 
this state is degenerate with the one that corresponds to the fourth entry. We therefore have to add the part of the 
perturbation that lifts this degeneracy: 



(n — l)a; 0 0 

0 (n — l)a; 0 

0 0 mo 



0 0 0 
0 0 0 


0 

0 



nuj 


0 0 
0 0 

0 0 
0 0 


0 {n + l)u; 0 

0 0 (n -I- l)u; 


(S12) 
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For large values of A this minimal addition to the Hamiltonian is indeed sufficient to obtain a good approximation 
for the quasienergies and Floquet states, because when A uj the Bessel functions are all much smaller than one, 
and we obtain two quasienergies with a separation given by the well-known expression AJq (2H/w). A problem arises, 
however, for small values of A/u). In the limit A/oj —)■ 0, one quasienergy is shifted up (from the unperturbed value) 
by A/2 and the other quasienergy is shifted down by the same amount. When that happens (and assuming that 
Lo and A are either exactly resonant or near resonance with each other), each one of the two quasienergies becomes 
degenerate or nearly degenerate with a similarly shifted level coming from the neighbouring pair of energy levels. As 
a result, if we now focus on the third entry (as above), we need to include the matrix elements in the Hamiltonian 
that describe mixing with the energy-level pair just above it in the matrix: 




(n — l)t 


A 


I v) 




f2A\ 


A , 

r2A\ 



h) 

0 



0 

0 


(n — l)o 


0 


A 




2A\ 

UJ ) 


2A 
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A (2A\ 




2A\ 
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I “) 

nu) 

0 
0 


(n -I- l)u; 0 
0 (n -I- l)u; 


V 


(S13) 
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In other words, after truncation and (with no loss of generality) choosing n = 0, we need to diagonalize the matrix 

( 


Hf. 


truncated — 



A , 

^2A\ 


A , 

(2A\ 

■UJ 

1 

to| 

V w / 

0 
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1 

K UJ ) 
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w J 

A /2A\ 
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2A\ 






2A\ 

w J 

_ , 2A 
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UJ 


0 

Aj f2A\ 
2 ^ \ U! J 
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We now perform a basis transformation H = HF,truncatedS, with 


5 = 


72 


/ 1 1 0 0 \ 
1-10 0 
0 0 11 
\ 0 0 1 -1 / 


and we obtain 


H = 


— 2uj — AJq 
0 
0 


f2A\ 



A /2A\ 


0 

0 



— 2uj -f AJq 
-AJi (— 1 


2A\ 

a 


, w / 

2A\ 


-AJi 


'2A" 

UJ 

a r I 2A^ 

— AJq ( - 

UJ 
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f2A\ 



A /2A\ 
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(S14) 


(S15) 


(S16) 


This matrix can be split into two 2x2 decoupled blocks, the relevant one of which is 


1 

2 


( -2uj + AJq (—\ 


UJ / 


A . /'2A 
—AJi — 


OJ 


—A(7i 
—AJo 



(S17) 
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The eigenvalues of this matrix (i.e. the quasienergies) are given by: 

U) 1 


2 2^ 


UJ — f — 


ei = 


ui 1 

-h 

2 2 


ui — A Jo ( — I 


1 2 


W yj 


A2 J2 ( — 

' UJ 


A2 J2 ( — 

' LO 


The Rabi frequency is given by the difference between the two quasienergies: 


Hr = 


^ r , 2A\ 

oj — A Jo ( — I 

w / 


A2J2 _ 


2^1 \ 


w / 


(S18) 


(S19) 


We note here that the expression for the quasienergy difference depends on our choice of quasienergies; choices other 
than the one used above would give expressions for the Rabi frequency that are different from the one given in Eq. (S19) 
by an integer multiple of w. The choice used here has the advantage that it gives the most natural expression for the 
Rabi frequency at small values of A, especially given the fact that there is only one dominant frequency component 
in the qubit population dynamics in the weak-driving regime. At higher values of A, it becomes more a matter of 
convention which frequency in the spectrum one defines as the Rabi frequency. In the weak-driving limit (A ^ w), 
this expression reduces to 


(w-A)VA2(^^^ , (S20) 

which upon assuming A « w reduces to the well-known expression 

Hr = ^(w-A)VA2. (S21) 

In the strong-driving limit {A » oj), Eq. (S19) reduces to 

Hr=uj- AJo , (S22) 

which, up to the physically unimportant differences of a shift by oj and a possible minus sign, is the well-known 
expression for the Rabi frequency in the strong-driving limit [S5, S6]. 

In Eig. S2, we compare quasienergies ij{A) at different driving frequencies obtained by numerically diagonalizing 
the Floquet Hamiltonian Eq. (S8) (solid lines) and the analytical formula Eq. (S18) (dotted lines). The numerical 
simulations use a Floquet Hamiltonian of matrix size of 101 x 101. The agreement between the numerical simulations 
and the analytical formula is very good in general. However, small disagreement is visible in the regime where A ~ w. 
This is because we ignore many off-diagonal matrix elements in Eq. (SIO), which contains Bessel function J„(2A/u;) 
with n > 2. We note that these Bessel functions vanish at A ^ w and A ^ w and they have maximum absolute 
values around A ~ w. 
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FIG. S2. Quasienergies versus the driving amplitude A for on-resonance driving condition, tu = A (a), and off-resonance driving 
conditions, cj = 1.4 x A (b) and cj = 0.6 x A (c). Solid lines are obtained by the numerical simulations, while dashed lines are 
obtained from the analytical formula. We show four neighboring quasienergies among the infinite number of solutions. The red 
and blue lines are quasienergies which correspond to the two inequivalent quasienergies, eq and ei respectively. 





